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ABSTRACT 

We present constraints on the non-linear coupling parameter /„; with the Wilkin- 
son Microwave Anisotropy Probe (WMAP) data. We use the method based on the 
spherical Mexican hat wavelet (SMHW) to measure the /„/ parameter for three of the 
most interesting shapes of primordial non-Gaussianity: local, equilateral and orthogo- 
nal. Our results indicate that this parameter is compatible with a Gaussian distribution 
within the two sigma confidence level (CL) for the three shapes and the results are 
consistent with the values presented by the WMAP team. We have included in our 
analysis the impact on /„/ due to contamination by unresolved point sources. The 



point sources add a positive contribution of A/, 



loc 

I 



2.5 ± 3.0, A = 37 ± 18 and 



A/°[* = 25 ± 14. As mentioned by the WMAP team, the contribution of the point 
sources to the orthogonal and equilateral form is expected to be larger than to the local 
one and thus it cannot be neglected in future constraints on these parameters. Taking 
into account this contamination, our best estimates for /„/ are —16.0 ^ f^f ^ 76.0, 
-382 s$ 202 and -394 < /°[* < 34 at 95% CL. The three shapes are com- 

patible with zero at 95% CL {2a). Our conclusion is that the WMAP 7- year data 
are consistent with Gaussian primordial fluctuations within 2a CL. We stress however 
the importance of taking into account the unresolved point sources in the measure- 
ment of /„/ in future works, especially when using more precise data sets such as the 
forthcoming Planck data. 
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1 INTRODUCTION 

During the period of inflationary expansion in the very early 
stages of the universe, primordial perturbations were gener- 
ated t hat are the seeds o f the structures that we can observe 
today ('Starobinskiri979': Guth'l98l':'Albrccht & Stcinhardt' 
[l982 : Lindc 1982. 1983.: Mukhanov ct al. 1992). These pri- 
mordial perturbations were linearly imprinted in the Cos- 
mic Microwave Background (CMB) anisotropies. Thus the 
study of the CMB anisotropies is a powerful way to under- 
stand the physics of the early universe. Many observational 
CMB projects, for example the NASA WMAlQ and ESA 
PlanclQ missions, dilTerent ground based 3D observational 
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campaigns of large scale structure and high energy accelera- 
tors are enabling us to understand better the properties and 
the evolution of the universe. From the several observational 
approaches that are available, the search for departures from 
Gaussianity in the CMB anisotropies with a primordial ori- 
gin has become a powerful way to discriminate among dif- 
ferent inflationary scenarios. Inflationary models such as 
the widely accepted standard, single-field, slow roll infiation 
predict low levels of non-Gaussianity whereas other models 
predict levels of non-Gaussianity that may be detected us- 
ing the data from current experiments iBartolo et al.l 20041 : 
lKomatsull2009l : lYadav fc Wandeltll2010l : lKomatsul l2010''). A 
detection of a deviation from Gaussianity with a primordial 
origin would rule out many inflationary models and would 
have far reaching implications in the physics of the early 
universe. 

The level of primordial non-Gaussianity is usually 
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parametrised by the non-linear coupling parameter /„; 
ll Verde et al.ll2000l : iKomatsu fc SpergeJ [2OO1I : iBartolo et all 
|2004 ). This parameter measures departures from zero in the 
values of the third order quantity known as the bispectrum, 
characterised through the shape function F{k\,k2,ks). 
The bispectrum is related to Bardeen's curvature per- 
turbations 'l?(k) through the 3-point correlation function 

($(ki)$(k2)$(k3)) = (27r)3<53(ki-Fk2-Fk3)F(fcl,fc2,fc3). 

Depending on the physical mechanisms of the different infla- 
tionary models the shape function can take different forms. 

In this paper we measure the levels of non-Gaussianity 
present in the WMAP data corresponding to the three par- 
ticular shapes (local, equilateral and orthogonal) that hav e 
been studied by the WMAP team ()Komatsu et al.ll201ll ). 
The shape function F(fci,fc2,fc3) of these types of non- 
Gaussianity, their CMB angular bispectra be-^e^e^ and the 
inflationary scenarios that generate these non-Gaussianity 
are described below. 

• Local shape. Significant non-Gaussianity of the local 
form can be generated for example in multi-field inflation- 
ary models (iKomatsu et al. 20051 : lKomatsull201Cf ). the cur- 



vaton model lifLvth et al.ll2003l ) ~he inhomo^eneou s reheat- 



mg scenario |Dvali et aI,P2004'; Bar tolo et al. [200^ '). models 
based on hybrid inflation (Lin 2009), etc. This sha pe is given 
by (see for example [C rcminclh ct al. 2006; Fcrgu sson et al.l 
l2010al : lYadav fc Wandelt.i201Q : iKomatsu et al.ii201l] ) 



F{ki,k2,k-i)=2A^Ui 



-1) 3_(„^_l) 



3-("=-l)r.3-(n,-l) 



+ 



,3-(n,-l), 3-{ne-l) 



(1) 



and its angular bispectr um is (see for example 
Fereusson erahl l2010a : .Yadav fc Wandehl l2010l : iKomatsui 



2OIOI ) 



X dx 



ae^{x)l3t^[x)l3e,^{x) + 
Pii {x)ae^ {x)l3e.^ (x) + jSi^ {x)l3e^ {x)ae^ {x 



(2) 



where A is the amplitude of the power spectrum P,s,{k) = 
Afe"'""*, Us is the spectral index and ae(x), Pe(x) 
are filter functions (see for example Komatsu fc Spergell 
2001 : IKomatsu et al ] |2005l : lFergusson et a ] |2010al : lKomatsul 



2OIOD . 

• Equilateral shape 

the equilateral form 



Significant non-Gaussianity of 
can be g enerated for example 



the Dirac-Bo r n-Infe l d infiation (ISilverstein fc Tone! 1200 



f-tolo ct ah 2004 
kani-Hamod ct al. 



Lan glois et al 
200i) 



2008 



ghost inflation 



sever al single-field infl ation- 

ary models in Einstein gravity llChen et al.l |2007|) etc 



This s hape is given by (see fo r example 



20061: Fereusson et al 



Komatsu et al. 



2OIII ) 



y (see to r example loremmet 
] l2010al : lYadav fc Wandeh 



Creminelli et al.l 
20ld: 



F{ki,k2,ks) = GA'Ui 
1 



3-(ns-l)l.3-(ns-l) 



3-(ns-l)^3-(n,-l) j.3- (n, - 1) ^3- (n, - 1) 



(fcifeA;3)2(4-".)/3 



.(4-n.) /3^2(4L.)/3^(4-n.) + P^^^^} 



fi,2 

and its angular b is pectrum is (see for example 



and Its angular b is pectrum is (see tor example 
Fergusson et al.l l2010al : lYadav fc Wandehl I2OI0I : IKomatsu! 



2OIOI ) 



dxx 



at^{x)l3i^{x)l3i^{x) + (2 perm) 



+I3i^(x)^i2{^)5i^{x) + (5 perm) - 25i^{x)5i^{x)5i^{x) 



(4) 



where 'ytjx) and 5i(x) are filter funct ions (see for example 
iFergusson eralll2010al : iKomatsullioiol ). 

• Orthogonal shape. Significant non-Gaussianity of 
the orthogo nal form can be generated in general single- 
field models (|Cheung et al.ll2008l : [Senatore et ahiboioh. This 



shape is given by (see for example Senatore et al. I2OI0I : 
lYadav fc Wandehil201Gl : IKomatsu et al.ll201lD 



F{kl,k2,k3) = 6AVni 

3 



3-(na~l)j,3-(n3-l) 



3-(n 



-1), 3-{n,-l) 



k 



+ {.(4-n. 



l.k-*-"a)/3> 2(4-n,)/3, (4-n,) 
ft,2 ft,3 



3-("s-l)i,3-(ns-l) 

8 

(fclfc2fc3)2(*-"=)/3 

-I- (5 perm) | 



(5) 



and its angular bi s pectrum is (see for example 



and Its angular bi s pectrum is ( 
lYadav fc Wandeltll201Gl : (Komatsullioiol ) 



dxx 



- (x)Pi>2 {x)pi^ {x) + (2 perm) 



8 , 



+I3t-^^{x)'yi^(x)5i^{x) + (5 perm) - -5i-^{x)5i^{x)5i^{x) 



(6) 



Many studies have been performed to constrain /„;, espe- 
cially for the local and the equilateral cases. The first con- 
straints on /„; were imposed using data sets with low resolu- 
tion or small sky coverage which led to large uncertainties in 
/„;. We can report anal yses using the Cosmic Backgr ound 
Explo rer (COBE) data (| Komatsu et a l. 2002; Cayon et al 



2003 ). MAXIMA data (|Cav6n et all 120031 : ISantos et al 



2003), the Very Small Ar ray (VSA) data (ISmith et al . 
20041 ) ■ the Archeops da ta (ICurto et al.| [2OO7I, I2OO8I) and 



the B OOMERang data (|De Troia et all l2007l : iNatoh et all 
[20l3). 

Once the WMAP data were available, significant im- 
provements were achieved in the precision of the estimation 
of fM 

Many studies have been developed to constrain the 
fni using WMAP data and based on different estimators. 
We can mention the different bispcc trum-based estima- 
tors (s ee for example Komatsu et all l2003l: Babich et al.l 



Creminelli et all 200( 



2004): IFergusson fc ShellardI l2007l: ISpergel et al 



Yadav fc Wandeh 
Komatsu et al. 



Fergusson & Shellard 

Smith etah 2009; Eisner fc Wandeh 2009 




— ■ Bucher et al 

2010! : iLiguori et ahilioiol : ISenatore et al.ir201Gl : ISmidt et al 



^ The improvement comes from a combination of large sky cover- 
age, high angular resolution and good sensitivity. This combina- 
tion improves the signal-to-noise ratio of /„; which for th e local 
case is proportional to log{tmax) llYadav fc Wandeltll201oh . 
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2OI0I: iFergusson et all l2010al lbl: iKomatsu et~all I2OIII : 

Fergusson fc ShellardI I2OIIII . The bispectrum is the most 



Gaussian a£_ coefficients: 



natural way to constrain /„; given its hnear dependence 
and the fact that in certain ideal conditions bispectrum- 
based estimators may be the optimal way to measure 
/„;. However, given that the data are contaminated by 
different non-Gaussian parasite signals and in most cases 
only a fraction of the sky can be used, it is convenient 
to use additional tools that can help to understand 
these effects better. We can mention the tests per- 
formed using the sp h erical Mexican hat wavelet (S MHW ) 
iMukheriee fc Wand 12004 ICurto et al.l l2009al l ^. l201iri . 



a HEALPix-based wavelet ( Casaponsa et alT 2011a b). 



a joint analysis with t he SMHW and neural networks 



Pietrobon et al. 20091: Rudiord et al. 20091: 


Pietrobon et al. 


201c 


; Rudiord et al.' 


20ldl:JPietrobon et al.ll2010l: Pietrobon 


201c 


: jCabcUa ct al. 


2010), the Minkowski functionals 



iMatsubaral I2OI0I: I Takeuchi et al.l 



2007I: Im 

20101). th 



__.|2010| ). the N-PDF distri- 

bution (IVielya fc Sanj|2009l. I2OI0I" or a B a yesian approach 



bution llVielya aanzllL^UUijl. I^Uimi or a Ij a yesian < 
l|Elsner et al.l l2010al : lElsner fc Wandeltl boiObh 



Other 

works use the 3D distribu tion of matter on large scales (see 



for example iDalal et al. 


2008: Matarrese fc Verdd 


2008; 


Slosar et al. 20081: Seliak 


1200a: .Desiacaues & Sel 


ak 


201c ; 


Xia et al.l 2010l: Baldauf et al.l 201ll: Hamaus et al. 


201 1|) to 



constrain the local /„;. 

In this paper we focus on the measurement of non- 
Gaussianity for the previous mentioned shapes using the 
estimator based on wa velets that has bee n f ormer l y use d 
JCurto et al.l l2009allbl. [ioTol . j201ll) ■ 
^2010al) 



to constrain local fni ^ .. . . 

We use the technique described bv IFergusson et al 
to produce non-Gaussian maps with the local, equilateral 
and orthogonal bispectra for WMAP resolution in realis- 
tic conditions of partial sky coverage and anisotropic noise. 
These maps are later used to evaluate the expected values 
of the wavelet third order moments aijk for each type of 
non-Gaussianity. We finally impose constrains on /„; for 
each shape using the wavelet estimator for the WMAP fore- 
ground reduced and raw data maps. As shown later, unre- 
solved point sources produce a significant bias in /„; that 
should be considered in the analyses of WMAP data and in 
the forthcoming analyses of Planck data, especially for the 
equilateral and orthogonal shapes. 

This paper is organized as follows. Section [2] presents 
the non-Gaussian maps that we have used to estimate the 
quantities needed for this analysis. In Section Owe present 
the method and the estimator used in this analysis to con- 
strain fni ■ The results of the analysis using WMAP data are 
presented in Section U and the conclusions are presented in 
Section (5] 



2 NON-GAUSSIAN SIMULATIONS 

Non-Gaussian Monte Carlo simulations are needed in or- 
der to calibrate the wavelet estimator. We have simulated 
ou r non-Gaussian maps fo llowing the algorithm described 
bv IFergusson et all (|2010al ). The non- Gaussian coeffi- 
cients can be written in terms of the bispectrum and the 



NG _ 1 \ ^ 

^2. "12,^3 ,"13 



£ £2 £3 
m m2 ma 



C«2 ^7^3 



(7) 



NG 



Using the fact that the shape functions of the local, equilat- 
eral and orthogonal bispectra are separable, we are able to 
reduce the number of sums in Eq. [T] This can be done in a 
straightforw ard way using Eg s. HI B]and|S]in Eq.[71 However , 
as stated by iHanson et al.l (|2009l ) ; Fergusson et al.l (|2010al ) , 
there are terms that may produce spurious divergences at 
low multipole s, large enough to affect the power spectrum of 
the final map. IFergusson et al.l l|2010ah located the divergent 
terms and provided equations for the local and equilateral 
shapes without these terms. A similar procedure can be per- 
formed with the orthogonal shape. In the next equations we 
present the non- Gaussian coefficients for each of the 

shapes without divergent terms. 

• Local bispectrum 

dxx^afXx) I d^nYi^{n)Mfi(x,n)Mp{x,n) (8) 
J 

Equilateral bispectrum 

dxx^^~3ae{x) J d^nYiln{n)Mi3{x,n)AIp{x,n) 

-2Se{x) J d^nYi^{n)Ms{x,n)Ms{x,n) 

+6'ye{x) J d^nYi^{n)Mi3{x,n)Ms{x,n)'^ 
Orthogonal bispectrum 

dxx^^~9ae{x) J d^nYl^{n)Mfi{x,n)Mp{x,n) 
~?<5i{x) j d^nYi^(n)Ms{x,n)Ms{x,n) 
-1-187^(2;) / d^ny/^(n)M^(a;,n)M^(x,n)| 



NG 



where 








ai{x) = 


2 

TT 


f°° 2 

/ k dkgTt{k)jt(kx) 
Jo 


(11) 


M^) = 


2 

TT 


poo 

/ k^dkP<t{k)gTe{k)je{kx) 
Jo 


(12) 


Jiix) = 


2 

TT 


/>oo 

/ edkPy\k)gTi{k)3t{kx) 
Jo 


(13) 


Si{x) = 


2 

TT 


poo 

/ edkPl'\k)gTi{k)ji{kx), 
Jo 


(14) 



greik) is the radiation transfer function that can be eval- 
uated using for example the CAMlfl or gTFasl[f| software, 
P#(fe) is the linear power spectrum, jiikx) is the spherical 



^ http://camb.info/ 

^ http:/ /gyudon. as. utexas.edu/~komatsu/CRL/indcx. html 
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Table 1. Quadrature in x integration used to compute a^J^\ 
We have used greater densit y of points near reioni z ation and re- 
combination as suggested by ISmith fc Zaldarriagal | |2006| ). Units 
for X are Mpc. 



the spherical harmonic space and by adding a Gaussian in- 
strum ental noise simulation in the real space (|Bennett et alj 
l2003h . 



^ a; ^ 9, 500 


64 points. 


Gauss-Legendre quadrature 


9,500 ^ a; ^ 11,000 


128 points 


, Gauss-Legendre quadrature 


11,000 < a; ^ 13,800 


64 points. 


Gauss-Legendre quadrature 


13, 800 a; ^ 14, 600 


170 points 


, Gauss-Legendre quadrature 


14,600 s£ a; ^ 16,000 


42 points. 


Gauss-Legendre quadrature 


16, 000 s£ X ^ 50, 000 


42 points. 


Gauss-Legendre quadrature 



Bessel function (of the first kind) and the M(x, n) maps are 
defined as 



Ma{x,h) = ae{x 
Mfi{x,h) ^^Pt{x) 



G yt,m(n) 

0-e,m 

yt,m{n) 



Ms{x,n) ^^3e{. 



, G y«,m(n) 



the 



power 

(aBB) 



(15) 
(16) 
(17) 

(18) 

spec- 



In Figure [T] we plot 

trum of the non-Gaussian terms _^ 

dxx^ae{x) J d'^nY(^{n)Mi3{x,n)Mfi{x,n), afj^'^'' = 
dxx^Si{x) J d^nYl^{n)Ms{x,n)Ms{x,n) and 

= dxx^je{x) J d'^nYi^{n)Mi3{x,n)Ms{x,n). 

The power spectrum of the Gaussian part is also plot- 
ted. We can see that these three terms add negligible 
extra-power to the full Gaussian plus non-Gaussian map. 
Once the terms are computed as a function of the 

bispectrum and the af^, the aim coefficients of a simulation 
with a given /„; can be written as airn ~ afm + fma'^m- 

In this paper we have generated a set of 300 non- 
Gaussian maps for the local, equilateral and orthogonal /„; . 
We have assumed a ACDM model using the parameter s 
that best fit the WMAP 7-year data (jKomatsu et al.ll201ll l. 
We have computed a power spectrum Ce and a transfer func- 
tion QTiik) using these parameters as inputs for the CAMB 
software (Lewis et al. 2000) up to Imax ~ 1535. The inte- 
grals in Eas.[8l[9land ll0l have been performed using a Gauss- 
Legendre quadrature. We have used a large density of points 
near reionization and recombination (see Table [T] for more 
details). A large number of points has been chosen in order 
to achieve convergence in the values of the Fisher matrix of 
the bispectrum. a% (fni ) for the three shapes . In Fi gure[2]the 
Fisher matrix ujp{f„i) jKomatsu &: Spergej|200ll 'l obtained 
with the three bispectra is plotted for different imax values. 
Note that these values are comparable w i th the values pre- 
sented for example bv lYadav Wandeltj (|2010l ). Once the 
atm of the simulations with non-Gaussianity are generated, 
we transform them into WMAP maps for each radiome- 
ter by convolving with the appropriate window functions in 



3 METHOD 

We use an estimator that is based on third-order statis- 
tics generated by the different possible combinations of the 
wavelet coefficient maps of the S MHW evaluated at certain 
angula r scales. See for example lAntoine fc VanderghevnsO 
J 19981): iMartmez-Gonzalez et all ( |2002l ): IVielval (|2007^ : 
iMartfnez-Gonzaled ( 20081 ) for detailed information about 
this wavelet. This estimator is described and used to search 
for blind non-Gauss ia n dev iations and constrain local /„; in 
ICurto et all (|2009al lbl. I2OIII ) . 

We consid e r the same set of angular scales Ri selected 
m ICurto et all l|201ll ). After evaluating the wavelet coeffi- 
cient map ■w{Ri;h) for each angular scale Ri we compute 
the third order moments qijk for each possible combination 
of thr ee angular scales {i, j, fc}. As mentioned in lCurto et al.l 
l|201lf ). the expected values of the cubic statistics are linearly 
proportional to fni 



(19) 



where the a^fc term is linearly related to the bispectrum. We 
evaluate these aijk quantities for the local, equilateral and 
orthogonal bispectra by averaging the values of the estima- 
tors obtained with the non-Gaussian simulations described 
in the previous Section. We then compute a statistic in 
order to constrain each /„; 



/ ohs 
(Qrst 



tfnl),m 



ijk,rst 

where q°j^ is the value of the statistics obtained for the ac- 
tual data map and C is the covariance matrix among the dif- 
ferent statistics (?ij4f|. The covariance matrix is estimated us- 
ing the Qijk statistics corresponding to 10,000 WMAP Gaus- 
sian s imulations. A detailed study described in lGurto et al.1 
l|201lD was carried out in order to compute correctly its 
inverse avoiding possible degeneracies. The Uijk statistics 
are estimated using the set of 300 non-Gaussian simulations 
transformed into WMAP V + W maps. Although we found 
analytical expressions f or the covariance matrix Cijk,rst and 
the aijk quantities (see lCurto et al.ll201l[ ). those expressions 
are only valid for the particular ideal case of full sky maps 
and white isotropic noise. For a realistic case, the analytical 
expressions become more complicated and the best practical 
approach to compute those quantities is using simulations. 

Finally, this estimator is also applied to a set of Gaus- 
sian maps in order to obtain an empirical estimate of the 
uncertainties of /„; . Additionally we also compute the value 
of the /„; Fisher matrix using the wavelet coefficients 



pi fni) 



Q-rst 



(21) 



° Note that for this estimator there is no need to subtract any 
Unear term due to the anisotropic noise as in the case of the 
KSW estimator. The reason is that the non-ideal aspects of the 
analysis (as the mask, the anisotropic noise, etc.) are included in 
the covariance matrix and the aijk coefficients. 
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power spectrum of 01^^" mops (local f^,— '^) 




power spectrum of o,^" mops (fni— "1 ) 



power spectrum of 0,^^^ maps (fm— 1 ) 




Figure 1. From left to right, the power spectrum of the non-Gaussian ^^f^^^^ ^^tm^^ coefRcicnts (lower line) compared 

with the Gaussian part of the power spectrum (upper hne). 



'^rCl'ni) fo'' local f^i 



Op(fni) for equilaterol f^. 




256 545 465 627 845 1139 1536 



738 887 1065 1279 1536 



""rCfni) orthogonal f^. 




512 615 738 BS7 1065 1279 1536 
multlpole I 



Figure 2. From left to right, the Fisher matrix ^^{f^-^i) versus Irnax for the local, equilateral and orthogonal bispectra b^^i^i^ defined 
in Eqs.[2][4]and[6] 



Table 2. Constraints on the local, equilateral and orthogonal f^i 
for the clean and raw maps and their uncertainties obtained with 
the Fisher matrix CTFifnl) with simulations (RMS). 



CASE 


raw /„; 


clean /„; 


0"F(/iii) 


MEAN 


RMS 


local 


25.0 


32.5 


22.5 


0.0 


23.00 


equilateral 


28.0 


-53.0 


145.0 


1.0 


156.0 


orthogonal 


-119.0 


-155.0 


106.0 


0.0 


112.00 



4 APPLICATION TO WMAP DATA 

4.1 Constraints on /„; using WMAP data 

We use the co mbined WMAP 7- y ear V and W band maps at 
the HEALPix (|G6rski et aLlbOOSl ) resolution oi N^e = 512. 
W e consider both raw a nd foreground reduced data maps 
as iKomatsu" et al.l (|201lh . The maximum multipole chosen 
in this analysis is 3A^sWe although the noise contamination 
starts to be significant at ^ ~ 1000. For the three shapes 
we find the best limits on /„; and provide the value of the 
Fisher and the simulated (j(f ni)- During all the analysis we 
use the WMAP KQ75 mask (|Gold et al.ll201ll ). In Tabled 
we summarize our results. We find that for the three cases, 
the parameters are compatible with zero at 95% CL. We 
would like to note the different effect that the foregrounds 
produce on different shapes: whereas it is negative for the 
local shape, it is positive for the equilateral and orthogonal 
shapes. For all the cases, arifni) is lower than the value 
obtained with simulations (~95% depending on the shape). 
We think that this small discrepancy is due to the limited 
number of simulations. We have checked that our estimator 
is unbiased. We have estimated the /„; values of 100 non- 
Gaussian simulations with an input /„; = 100 and used 
the remaining 200 non-Gaussian simulations to estimate the 



Oijk- The results are /^f = 99.5 ± 29.5, f^l = 98 ± 150 and 
/nF* = 97 ± 118, which are clearly compatible with the input 
fni taking into account the expected errors in the mean for 
the available number of realizations. Our best estimates for 
the clean maps are: 

• LocaO: /„; = 32.5 ± 22.5 (68% CL) 

• Equilateral: /„; = -53 ± 145 (68% CL) 

• Orthogonal: f„i = -155 ± 106 (68% CL) 



results presented by 
one sigma error-bars. 



The values m atch well the 
iKomatsu et al.l (|201ll ') within 
The differences can be explained by the different sensitivity 
of the bispectrum and wavelet estimators to the possible 
non-cosmological residuals present in the data. 



4.2 Point source contribution 

We have also estimated the contribution of undetected point 
sources using the sou rce number counts dN/dS derived from 
Ide Zotti et all l|2005h . We have used point source simula- 
tions based on this dN/dS. We have chosen a maximum 
flux for the bright sources such that the power spectrum 
for the Q band is compatible with the value provided by 



the WMAP team, 



-^7)3 — 



0.0090 ± 0.0007 iiK^ sr in an- 



tenna units ()Larson et al.l [loij) . We have estimated the 



^ Using a set of 30 non-Gauss ian simulations generated following 
the procedure by iLiguori et al. (2003, 2007) our best estimate 
is fnl = 37 ± 21 (68% CL). These maps were generated in a 
different way: the non-Gaussianity is introduced in the primordial 
curvature perturbation "I>(x) = ^^(x) -|- /„; (<3>|^(x) — ($|^(x))^ 
and then extrapolated to the CMB. This process add extra non- 
Gaussianity at higher moment s where as the procedure used in this 
paper and in lFergusson et al.l (|201o3) just adds non-Gaussianity 
to the third order moments (bispectrum). 
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best-fitting /„; value for two sets of 1,000 maps. Tlie first set 
consists of 1,000 Gaussian CMB + noise maps and the sec- 
ond consists of the same Gaussian CMB -I- noise maps plus 
the point source maps. For each map with point sources 
we estimate its best-fitting /„; parameter and compare with 
the value obtained for the same map without point sources. 
The difi'erence A/„i provides an estimate of the impact on 
fni due to the unresolved point sources. The point sources 
add a contribution of A/„; = 2.5 ± 3, A/„; = 37 ± 18 and 
A/„i = 25 ± 14 for the local, equilateral and orthogonal 
forms respectively. 

To check further these results, we have used an alterna- 
tive method to estimate the point source contamination to 
fni given by the expression 

^Ui = -F^ p;ri , (22) 

where {qijk)pa is the expected value of the third order mo- 
ments due to the point sources. The results are A/„; = 2.5, 
A/ni = 38 and Af^i — 24 which agree with the values pre- 
viously obtained with simulations. Taking into account the 
point source contribution, our best estimates of /„; are: 

• Local: /„; = 30.0 ± 22.5 (68% CL) 

• Equilateral: f„i = -90 ± 146 (68% CL) 

• Orthogonal: /„; = -180 ± 107 (68% CL) 

Fig. |3] contains the histograms of the best-fitting /„; values 
for each shape corresponding to 1,000 CMB -I- noise Gaus- 
sian simulations and the values of the data after the point 
source correction. Note that the point sources add a signifi- 
cant contribu tion to the equilateral and orthogonal shapes. 
We agree with lKomatsu et all (|201ll ) that the WMAP seven- 
year data are consistent with Gaussian primordial fluctua- 
tions for the three considered shapes. Planck will be able 
to address this issue with more detail due to its increased 
sensitivity and power to clean the signal. 



5 CONCLUSIONS 

We have imposed constraints on primordial non-Gaussianity 
with the WMAP 7-year data using the wavelet based esti- 
mator. In this analysis we have considered the combined 
V-l-W maps and the KQ75 mask. In particular, we have fo- 
cused in three shapes with particular interest for the physics 
of inflation in the early universe: the local, equilateral and 
orthogonal bispectra. 

We have simulated the non-Gaussian maps for each 
of the considered shapes and estimated with these simu- 
lations the required quantities for our estimator. Our re- 
sults are compatible with the values obtained by the WMAP 
team and our uncertainties are very similar to the er- 
ror bars obt ained with the optimal bispectrum estimator 
(jKomatsu et al. 2011). 

In addition we have estimated the contribution of the 
point sources. In the particular case of the local fni , the con- 
tribution is A/^""^ = 2.5 ± 3.0. T his is similar to the v alues 
obtained by the WMAP team l|Komatsu et al.l [2OI1I ) and 
its contribution to the parameter is not signiflcant. How- 
ever we have detected a non-negligible contribution to the 
equilateral and orthogonal shapes due to the unresolved 
point sources. In particular, we have found A/^^ = 37 ± 18 



and A/°r* = 25 ± 1 4 (68%CL). These la rge values were 
already predicted by iKomatsu et al.l (|201ll ) although they 
did not provide actual flgures. This contribution should be 
taken carefully into account in future constraints on f„i with 
WMAP and Planck data. Considering the point sources, our 
best estimates of /„; are flT = 30.0 ±22.5, = -90 ±146 
and /°p = —180 ± 107. The three shapes are compatible 
with zero at 95% CL. Our conclusion is that the fni param- 
eters are compatible with ze ro within the 2 a CL and our 
results are in agreement with iKomatsu et all (I2OIII ). 

The wavelet estimator has been tested and carefully 
checked with th e available W MAP d ata in this and several 
previous works (jCurto et al. 2009a b, 1201 ll ). It is now ready 
and being upgraded to analyse the forthcoming Planck data. 
In future works we will also use the wavelet estimator jointly 
with ne ural networks to constrain these shapes along the 
lines of ICasaponsa et al.l l|2011bl ) where this procedure has 
been already applied for the local /„; using WMAP data. 
This later process helps to speed up the calculations since 
it is not necessary to estimate the covariance matrix of the 
cubic statistics and it avoids all the possible complications 
in the computation of the inverse covariance matrix. 
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